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We investigate the possibility that dark energy is scaling with epochs. A phenomenological model 
is introduced whose energy density depends on the redshift z in such a way that a smooth transition 
among the three dominant phases of the universe evolution (radiation era, matter domination, 
asymptotic de Sitter state) is achieved. We use the WMAP cosmic microwave background data 
and the luminosity distances of Type la Supernovae to test whether the model is in agreement with 
astrophysical observations. 
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INTRODUCTION 

One of the greatest challenge in modern cosmology is to 
identify the nature of the dark energy component which 
is causing the observed accelerated expansion of the uni- 
verse. Since a cosmological constant, while in agreement 
with current observations, is theoretically flawed, several 
alternatives have been proposed. Some of the popular 
candidates to explain the observations are a slowly-rolling 
scalar field, "quintessence" or a "k-essence" scalar 

field with non-canonical kinetic terms in the Lagrangian 
[Sj-j^J or "coupled quintessence" where the scalar field is 
non-minimally coupled with gravity 3] . The cosmologi- 
cal acceleration can be also achieved considering geomet- 
rical counterparts in the gravitational Lagrangian other 
than the standard Ricci scalar of General Relativity. This 
fact allows to define a sort of effective curvature pressure 
and curvature energy density which act as a time- varying 
cosmological constant . 

One way to distinguish whether the dark energy is due 
to a cosmological constant, to a scalar field or something 
else is to measure the equation of state, wx, the ratio 
of the pressure px to the energy density px ■ A cosmo- 
logical constant always has wx = — 1 while scalar fields 
or curvature counterparts generally have an equation of 
state which differs from unity and varies with red shift z. 
Through measurements of type la supernovae (SN-Ia), 
large-scale structure and Cosmic Microwave Background 
(CMB) anisotropics, the equation of state may be deter- 
mined accurately enough in the next few years to find 
out whether dark energy is actually different from a cos- 
mological constant or not. 

An important general property of these so-called 
'tracker' models (and of any dark energy model which 
aims to alleviate the fine tuning problem of the cosmo- 
logical constant) is that the scalar field equation of state 
(and its energy density) remains close to that of the dom- 
inant background component during most of the cosmo- 
logical evolution. 



For example, in power-law potential like V = 
the equation of state generally remains closer to the back- 
ground value Wx = awBl{ci + 2) while the ratio of the 
energy density of the scalar field to that of the domi- 
nant component gradually increases. In models based on 
exponential potential V = K)e^^"^ (0.11) wx(z) mim- 
ics exactly the scaling of the dominant background in 
the attractor regime {wx = wb) and if the background 
component scales as pB — Po(^)"; then the scaling field 
approaches an attractor solution, and its fractional en- 
ergy density is given by U,x = constant with redshift. 
As further example, in "k-essence" models, the fc-essence 
undergoes two transitions in its behavior, one at the on- 
set of matter-domination and a second when fc-essence 
begins to dominate over the matter density. During the 
radiation-dominated era, the fc-essence energy tracks the 
radiation, falling as l/a^ where a is the scale factor. The 
onset of the matter-dominated era automatically triggers 
a change in the behavior of fc-essence such that it begins 
to act as an energy component with wxiz) < 0. When 
fc-essence overtakes the matter density, wx{z) changes 
to another value around —1, the precise value of which 
depends on the detailed model. 

Given the new release of cosmological data from high- 
precision measurements of CMB aniso trop ies (see e.g. 
'23'|) and SN-Ia luminosity distances (|25j) (which are 
now providing a ~ 18a evidence for a dark energy com- 
ponent) is therefore extremely timely to check if any hint 
for a "scaZmi?" dark energy is present in the data. 

Moreover, recent analysis of SN-Ia data (see e.g. 
and with model independent parameterization have 
found that dark energy which evolves with time provides 
a better fit to the SN-Ia data than a standard cosmolog- 
ical constant. 

In this paper we use a phenomenological approach to 
constrain a dark energy component with an evolution- 
ary behavior similar to the models mentioned above. In 
particular, we use a toy model whose energy density de- 
pends on the redshift z in such a way that a smooth 
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transition among the main three cosmological scahng 
regimes (radiation, matter and dark energy) is achieved 
M- We then use the Wilkinson Microwave Anisotropy 
Probe (WMAP) CMB data and the luminosity distances 
of SN-Ia to test whether the model is in agreement with 
astrophysical observation and/or any evidence for "scal- 
ing" dark energy is present in the data. 

This approach has the main advantage of being the- 
oretically well motivated since a scaling model can rea- 
sonably approximate the behavior of most of the dark 
energy theories on the market (cosmological constant in- 
cluded). An analysis with an higher number of parame- 
ters to describe the dark energy evolution would probably 
be better suited to detect variations from a cosmological 
constant or to test models with a rapidly evolving equa- 
tion of state. However, allowing more degrees of freedom 
could introduce serious degeneracies, fit unknown sys- 
tematics and produce final results of difiicult theoretical 
interpretation. The analysis and results presented here 
can be therefore considered as complementary to recent 
analysis which sampled a wider set of parameters (see e.g. 
^3]) or on the contrary restricted the study to a constant 
with redshift equation of state (see e.g. |2^, 0|). The 
plan of the paper is as follows: in the next section, the 
phenomenological model is illustrated in detail. In Sec- 
tion 3, the analysis method is briefly discussed while in 
Section 4 we present the results. Section 5 is devoted to 
conclusions. 



A PHENOMENOLOGICAL MODEL FOR 
SCALING DARK ENERGY 

Let us now illustrate in detail the phenomenological 
method used for our analysis. The energy density of our 
scaling dark energy model evolves with redshift as (see 

My- 



Px{z) = A 1 + 
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where A is a normalization constant, related to the to- 
day dark energy density Vlx, and Zs and z^ are free 2 
parameters that identify the three epochs of scaling. The 
corresponding dark energy equation of state is indeed: 
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which depends on the parameters Zg and such that: 



Wx 



1/3 for z» Zs 



Wx ~ for Zb << z << Zs , 

Wx ^ —1 for z << Zb . 

Therefore the model we obtain is able to mimic a fluid 
following flrst a radiation equation of state, then a mat- 
ter phase and flnally approaching a deSitter phase with 
constant energy. This model is useful to identify and/or 
constrain a cosmological imprint of scaling dark energy 
in the data. The same model can be further extended to 
a fluid with transitions between more generic equations 
of state. 

In Fig. n] the behavior of the energy density of the 
3 components in our model (matter, radiation and dark 
energy) in function of the scale factor are plotted. As 
one can see, if the matter to radiation transition redshift 
is much smaller than the A-CDM redshift of equivalence 
Zs < ZEq ^ 3200, the dark energy can be again domi- 
nant in the past. However this can be tuned by Zb which 
shifts the dark energy transition between a cosmological 
constant and matter. For example, a radiation to matter 
transition at redshift Zg ^ 100 can still be in agreement 
with a negligible dark energy contribution in the past, 
providing that Zb > 5. 

The time variation in wx is small in comparison to 
the expansion rate of the universe. We assume purely 
adiabatic contributions to the perturbations in the spec- 
trum. The sound speed is therefore fully determined by 
wxiz) and we integrate the evolution equations for the 
density and velocity (we neglect shear) perturbations in 
the dark energy fluid as in While the adiabatic 

approximation is not completely correct for most scalar 
field or quintessence models, since there dissipative pro- 
cesses generate entropic perturbations in the fluid and a 
more general relation is needed, such a model can be con- 
sidered anyway as a good start for analyzing the current 
data. 

In Fig. 121 top panel, we plot several power spectra com- 
puted with CMBFAST ^18] in function of Zb with mat- 
ter density Qm = 0.35, f2x = 0.65 and the dark energy 
matter-radiation transition redshift fixed at Zs = 5000, 
well after the redshift of equivalence in standard cold 
dark matter model dominated by a cosmological constant 
(A-CDM). As we can see, increasing has the effect of 
mimicking more and more the cosmological constant be- 
havior. On the other hand, decreasing Zb increases the 
effective equation of state Wef / shifting the peaks in the 
CMB spectrum toward smaller angular scales (see e.g. 

In Fig. 121 bottom panel, z;, = 5 is fixed and the depen- 
dence of Zs is studied. In this case, the effect is smaller, 
and it is clear that as soon as Zg becomes smaller than 
the redshift of equivalence in A-CDM, fix can dominate 
again as a relativistic component and leave an imprint on 
the CMB spectrum through the Early Integrated Sachs- 
Wolfe effect (see e.g. [l3|). 
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FIG. 2: Top Panel: CMB power spectra normalized at £ = 111 
for a universe with flm = 0.35 and fix = 0.65. The transition 
redshift matter-radiation has been fixed to Zs — 5000. Models 
with different values of zi, are plotted: red {zt = 1), green 
{zt — 2), blue {zi, — 5), light blue {zt — 20). Bottom Panel: 
CMB power spectra normalized at £ = 111 for a universe 
with Qm ~ 0.35 and fix ~ 0.65. The matter-dark energy 
transition redshift has been fixed to zt — 5. Models with 
different values of Zs are plotted: green {zs — 50), blue {zs = 
100), yellow {zs = 500), red {zs = 5000). The WMAP data 
points are also plotted for comparison. 



FIG. 1: Evolution of the overall energy density of the Universe 
with redshift in difi'erent theoretical frameworks. In the pan- 
els, the phenomenological dark energy model described in the 
text has Zs — 100 (grey line) and Zs = 1000 (black line). The 
top, middle and bottom panels show the cases with Zb = 1, 
Zb = 2 and Zb = 5 respectively. In all plots, the short dashed 
line is the matter component while the long-dashed line is 
radiation (photons and 3 massless neutrinos) contribution. 



THE ANALYSIS METHOD 

In order to bound this phenomenological dark energy 
model, we consider a template of flat, adiabatic, X-CDM 
models computed with CMBFAST }^ . We sample the 
relevant parameters as follows: ^cdmh'^ = 0.05, ...0.20, 
in steps of 0.01; ^bh? = 0.015, 0.030 (motivated by 
Big Bang Nucleosynthesis), in steps of 0.001 and ilx = 
0.05, ...,0.95, in steps of 0.02. 

We sample the 2 parameters of the dark energy model 
in the range Zb = 0.05, ...,7.55, in steps of 0.5 and Zg = 
20,. ..,740 in steps of 80. 



The value of the Hubble constant in our database is not 
an independent parameter, since it is determined through 
the flatness condition. The conservative top-hat bound 
0.55 < h < 0.85 is adopted and the la constraint on 
the Hubble parameter, h = 0.71 ± 0.07, obtained from 
Hubble Space Telescope (HST) measurements 0, is also 
considered. 

We allow for a reionization of the intergalactic medium 
by varying the Compton optical depth parameter Tc over 
the range Tc = 0.05, 0.25 in steps of 0.02. 
For the CMB data, we take into account the recent tem- 
perature and cross polarization results from the WMAP 
satellite ('is']) using the method explained in ([l^l) and 
the publicly available code on the LAMBDA web site. 
Finally constraints obtained from the luminosity mea- 
surements of Type la Supernovae (SN-Ia) are incorpo- 
rated from Hij using the GOLD dataset. The SN-Ia 
luminosity data are independent of the neutrino energy 
density but, as the CMB data, are helpful in breaking de- 
generacies between the parameters we are going to con- 
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sider. 

THE RESULTS 

In Figure O the likelihood contours are plotted in the 
{zs, Zft) plane only for the CMB data. As we can see 
the current data doesn't show any evidence for scaling 
dark energy and we can only derive weak lower limits on 
the two parameters of the model. We have found that 
Zs > 60 and zt > 3.8 at 1 — a. From the same Figure is 
clear that an interesting correlation exists between the 2 
parameters. Namely, the current CMB data do not fa- 
vor or provide evidence for an extra matter or radiation 
component: if Zg is too small then it is necessary to con- 
sider larger values of Zb in order to have the extra dark 
energy component not dominant in the past. At high z^, 
we have an asymptotic value 95% c.l. of zi, > 3.0. Also 
plotted on the figure are the likelihood contours derived 
by a combined CMB-|-SN-Ia analysis. The inclusion of 
the SN-Ia data improves the constraints to Z(, > 5 and 
Zs > 100 at 1 — (7. SN-Ia data are insensitive to variations 
in Zs since they are probing only redshifts z < 1.5, but 
provide anyway complementary constraints on Zb and the 
matter density Qm- 

This is better explained in Figure 01 where we over- 
impose the likelihood contours in the (fim, Zb) plane 
from CMB and SN-Ia analysis. As we can see, the cur- 
rent SN-Ia data does not provide evidence for dark en- 
ergy evolution. However, low values of Zb are compatible 
with the SN-Ia data if one decreases the amount of the 
matter component. The SN-Ia data is consistent with 
rim = but in this case the dark energy model behaves 
like an unified dark energy model with transition red- 
shift Zb ^ 0.5. On the other hand, as we can see from the 
plot, lower values of Zb are compatible with CMB data if 
one increases the matter density. This is easily explained 
from the fact that a lower Zb results in higher values for 
the effective dark energy equation of state. The direction 
of degeneracy in the plot in the case of the CMB data 
is therefore only a consequence of the geometrical degen- 
eracy present in angular diameter distance data at high 
redshift. The lower limit on Zb from the CMB data comes 
mainly from our assumptions on the possible values of 
the Hubble parameter. Combining the CMB and SN-Ia 
further breaks this degeneracy, improving the lower limit 
on Zb and excluding at high significance an unified dark 
energy model with — 0. 

CONCLUSIONS 

A phenomenological scaling dark energy model is dis- 
cussed in this paper matching it with the current CMB 
and SN-Ia data to identify the signatures of a possible 
cosmological evolution of dark energy, as expected and 
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FIG. 3: Likelihood probability contours at 68% (light blue), 
95% (blue) and 99% (dark blue) in the Zs-Zb plane from 
WMAP. The solid,long dashed and short-dashed lines are the 
68%, 95% and 99% from a WMAP-f-SN-Ia analysis. 
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FIG. 4: 68% and 95% likelihood probability contours in the 
Qjn-zt plane from WMAP and SN-Ia. 

predicted in several theories. We do not take into account 
particular scalar field or quintessence models but discuss 
how the state equation of cosmic fluid, depending on red- 
shift, scales with respect to the epoch passing from a ra- 
diation regime to a dark-energy-matter-dominated era. 
The approach is extremely general since the dynamical 
behavior which we discuss should be the one expected for 
the most of the models present in literature. 

We found that the current data does not show evidence 
for cosmological evolution of dark energy, providing the 
68%c.l. bounds Zb > 5.0 and Zs > 100 constraining the 
presence of scaling dark energy in the universe. A simple 
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but theoretically flawed cosmological constant still pro- 
vides a good fit to the data (see also the discussion in 

i). 

Interesting correlations between the parameters used 
in the analysis of single datasets are present. For example, 
lower values of the matter density (ri,„ ~ 0.1) could hint 
for an evolution in the SN-Ia data, but these values are 
ruled out by CMB measurements. 

Dark energy models with a subdominant contribution 
to the overall energy density of the universe for most of 
the cosmological evolution are clearly preferred. While 
this condition is easy to achieve for many models based 
on scalar fields or topological defects (see e.g. ^^), the 
quantity of information we can hope to extract from fu- 
ture data about dark energy is more limited. 

However, our approach could be further improved 
and extended to more general scaling solutions. In 
this respect, the incoming wide cosmological surveys as 
PLANCK or SNAP will provide essential data. 
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